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We discuss theoretical models for the cooperative binding dynamics of ligands to substrates, such 
as dimeric motor proteins to microtubules or more extended macromolecules like tropomyosin to 
actin filaments. We study the effects of steric constraints, size of ligands, binding rates and inter- 
action between neighboring proteins on the binding dynamics and binding stoichiometry. Starting 
from an empty lattice the binding dynamics goes, quite generally, through several stages. The first 
stage represents fast initial binding closely resembling the physics of random sequential adsorption 
processes. Typically this initial process leaves the system in a metastable locked state with many 
small gaps between blocks of bound molecules. In a second stage the gaps annihilate slowly as 
the ligands detach and reattach. This results in an algebraic decay of the gap concentration and 
interesting scaling behavior. Upon identifying the gaps with particles we show that the dynamics 
in this regime can be explained by mapping it onto various reaction-diffusion models. The final ap- 
proach to equilibrium shows some interesting dynamic scaling properties. We also discuss the effect 
of cooperativity on the equilibrium stoichiometry, and their consequences for the interpretation of 
biochemical and image reconstruction results. 

PACS numbers: 68.45Da, 82.20Mj, 87.16Nn 



I. INTRODUCTION 

How does a system evolve towards its steady state? 
Sometimes the answer is quite simple and the relaxation 
process is merely an exponential decay. If the deviations 
from equilibrium are small Onsager's regression hypothe- 
sis Q asserts that the relaxation is governed by the same 
laws as the fluctuations in equilibrium. This hypothe- 
sis certainly fails for systems with an absorbing steady 
state such as simple models for diffusion-limited chemi- 
cal reactions [|, |], [|, [|, g|. Here there are no fluctua- 
tions in the steady state but the approach towards the 
absorbing state is critical in the sense that it exhibits 
slow power-law decay and universal scaling behavior 0. 
Most of these models are chosen to be mathematically 
transparent hoping that they still resemble some of the 
essential features of actual systems occurring in nature. 
Unfortunately, experimentally accessible systems where 
the above theoretical ideas can be tested explicitly have 
remained rare to date. 

In this contribution we discuss the kinetics of some 
macromolecular assembly processes relevant for the for- 
mation of functional structures in cells. In particular, we 
are interested in the dynamics of ligand-substrate bind- 
ing, where the substrate is a one- or two-dimensional lat- 
tice and the ligands are dimers or oligomers. Examples 
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for such systems are illustrated in Fig. [lj. Binding of 
dimeric myosin to actin filaments j|] can behave in a 
similar way as dimeric kinesin on microtubules. Simi- 
lar interactions of large supramolecular biological poly- 
mers with protein ligands such as DNA with proteins or 
viruses with antibodies are also a quite intensive area of 
research. As will become clear in later chapters the kinet- 
ics of these systems shows interesting anomalous dynam- 
ics which is closely related to the mathematical models 
of simple diffusion-reaction systems discussed above. 

Fig. |l|a shows a schematic representation of a "decora- 
tion experiment" where dimeric motor enzymes (ligands) 
are deposited on their corresponding molecular track 
(substrate) |], [n|. For kinesin motors, these tracks are 
microtubules, hollow cylinders usually consisting of 13 
protofilaments, linear polymers composed of alternating 
a- and /3-tubulin subunits. The kinesin binding sites are 
located on the (3 subunits (dark spheres) which form a he- 
lical (wound-up rhombic) lattice with a longitudinal pe- 
riodicity of 8 nm. Kinesin is a mechanochemical enzyme, 
which transforms (through an isothermal stochastic pro- 
cess) the chemical energy obtained from the hydrolysis 
of adenosine triphosphate (ATP) into motion along mi- 
crotubules. Roughly speaking, kinesin has the follow- 
ing building plan: two globular heads (also called motor 
domains) with the functionality of "legs" are joined to- 
gether by a coiled-coil region into a tail which can bind 
to some cargo. The typical size of these proteins is in the 
range of several tens of nanometers. 

Decoration techniques are usually performed in the ab- 
sence of ATP hydrolysis. Then the motor enzymes are 
"passive" ligands which bind and unbind from their tubu- 
lin binding sites but do not actively move along micro- 
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FIG. 1: a) Kinesin dimers can bind one head or both heads to a one-dimensional lattice (tubulin protofilament). Binding 
sites are located on /3-tubulin subunits (dark) while a-subunits (bright) are irrelevant for our study, b) Tropomyosin binds to 
an actin filament by occupying 7 lattice sites. There is a strong attractive interaction between the ends of bound tropomyosin 
molecules. j— i 

b) © 2001, Biophysical Society jl4 



tubules. These systems have traditionally been used in 
biophysical chemistry to investigate the structure and the 
binding properties of kinesin ^J, i.e. after waiting 
for the system to equilibrate the binding stoichiometry 
is measured and the structure is determined by cryo- 
electron microscopy followed by 3D image reconstruc- 
tion. These investigations are key elements in under- 
standing the binding patterns of kinesin motor domains 
under changing nucleotide conditions to get a complete 
picture on the different conformational states, which are 
involved in the kinesin walking mechanism flllfl . 

In the following sections we present a theoretical anal- 
ysis of dimer adsorption-desorption kinetics with com- 
peting single and double bound dimers. Such analysis is 
extremely important for a quantitative analysis of deco- 
ration data fl^j ; a theoretical analysis gives the binding 
stoichiometry in the equilibrium state in terms of bind- 
ing constants for the first and second head of the dimer 
molecule. The dynamics of the approach to equilibrium 
is useful to estimate when an experimental system can be 
considered as equilibrated. Even more importantly, time- 
resolved decoration experiments (e.g. by using motor en- 
zymes labelled by some fluorescent marker) combined 
with our theoretical analysis could provide new informa- 
tion about reaction rates which are to date not known 
completely. Understanding the kinetics of passive mo- 
tors is undoubtedly a necessary prerequisite for studying 
the more complicated case of active motors at high den- 
sities The model is also interesting in its own right 
since it contains some novel features of non-equilibrium 
dynamics of dimer adsorption-desorption models. 

Fig. |l|b shows a schematic representation of 
tropomyosin binding to an actin filament Actin 
filaments are one of the major compontents of the 
cytoskelcton. Like microtubules, they contribute to the 
mechanical stability of the cell and serve as tracks for 



molecular motors from the myosin family. This function 
is especially pronounced in muscle cells where both 
actin and myosin form filaments that can slide between 
each other and thereby cause the muscle to contract. 
Tropomyosin plays the key role in the control of skeletal 
muscle contraction. It binds to actin along its binding 
sites for myosin motors. When calcium ions are released 
as a response to a nerve signal, they cause tropomyosin 
to shift laterally thereby clearing the binding sites and 
allowing myosin to bind and produce force. Tropomyosin 
binding to action has several features which make it 
an interesting model system of statistical physics; each 
molecule covers seven actin monomers and interacts 
strongly with other tropomyosin molecules. As a 
consequence gaps between bound molecules can take 
a long time to heal. Although the gap dynamics has 
not yet been measured experimentally, the amount of 
data gathered in other kinetic studies provides plenty 
of information on the model parameters and allows to 
make quantitative predictions about the relaxation time. 
As it turns out ]l4| ], these relaxation times can be as 
large as hours, which is short enough compared to the 
lifetime of actin filaments in a muscle cell, but essential 
when planning experiments with actin and tropomyosin 
assembled in vitro, especially when studying the myosin 
regulation |]l5| . 

The outline of this article is as follows. In Section || 
we define the model for decoration of microtubules with 
dimeric motors which can bind either with one or with 
two heads as first introduced in Ref. [jl2|. In Section III 
we determine analytically the equilibrium state of this 
model. We study the dynamics of the model in Sec- 
tion |v|, which extends the results presented in Ref. jl6) . 
Section ^ studies the dynamics of the two-dimensional 
model with interactions and Section [v| the dynamics of 
the one-dimensional fc-mer model, which is a more gen- 
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eralized version of the results described in Ref . H] . In 
the final section we give a summary and an outlook on 
future challenges in the field. 



II. DEFINITION OF THE MICROTUBULE 
DECORATION MODEL 

Our model describes the experimental situation as 
found in most decoration assays. It starts with an empty 
tubulin sheet surrounded by a solution of double-headed 
kinesin molecules. The kinesin dimers can either attach 
with one or two heads onto binding sites located on (3- 
tubulin. There seems to be convincing evidence that ki- 
nesin heads can bind on two adjacent binding sites only 
in longitudinal but not in lateral direction 
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introduces a strong uniaxial anisotropy and distinguishes 
the adsorption process of protein dimers from simple in- 
organic dimers. The attached heads can also detach at 
some rate. A schematic representation of the decoration 
process onto tubulin sheets is given in Fig. |2[ 




FIG. 2: Schematic representation of the binding of kinesin 
dimers to a tubulin sheet. The binding sites are located on the 
/3-tubulin subunits marked as grey balls. A dimer can attach 
either with one head to one binding site or with two heads 
on two adjacent sites along the same protofilament. Each 
binding site can be occupied by at most one kinesin head. 
Neglecting interactions between kinesin molecules results in 
a one-dimensional model with dimers decorating individual 
protofilaments. 

To begin with, let us take into account only steric 
(hard-core) interactions and for now neglect nearest 
neighbour attractive interactions. Then we are left with 
a one-dimensional problem of kinesin dimers decorating 
a single protofilament (one-dimensional lattice) as shown 
in Fig. || and defined as follows. p9| Kinesin is considered 
as a dimeric structure with its two heads tethered to- 
gether by some flexible joint. Hence each dimer (kinesin 
protein) can bind one of its two heads (motor domains) 



to an empty lattice site |T|]. The binding rate fc+i c for 
this process is proportional to the solution concentration 
c of the dimeric proteins. Successively, the dimer may ei- 
ther dissociate from the protofilament with a rate or 
also bind its second head to an unoccupied site in front 
(f) of or behind (b) the already bound head (with front 
we always refer to the direction pointing towards the "+" 
end of the microtubule, which is the direction of motion 
for most motors from the kinesin family). Since kinesin 
heads and microtubules are both asymmetric structures 
the corresponding binding rates k +2 an d k b +2 are in gen- 
eral different from each other. The reverse process of 
detaching a front or rear head occurs at rates k_ 2 and 
k_2- A reaction scheme with all possible processes and 
their corresponding rate constants is shown in Fig. |i| In 
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FIG. 3: Reaction scheme for all binding and unbinding pro- 
cesses. The rates for binding and unbinding of the first head 
are given by fe+i c and k-i, respectively. The second head 
may bind to an unoccupied site in front (f) of or behind (b) 
the already bound head with rates k^_ 2 and k h +2 - The cor- 
responding reverse process of detaching a front or rear head 
occurs at rates fc^ 2r and k b _2- 

© 2001, Academic Press |l2j 

typical decoration experiments there is no external en- 
ergy source, i.e. no ATP hydrolysis. Then the binding 
rates are not all independent of each other, but detailed 
balance dictates that the ratio of on- and off-rates has to 
equal the equilibrium binding constants 



k+i 
k-i 



and K 2 



k b 

"■+2 



"+2 



(1) 



A particular coverage of the lattice is described as a se- 
quence of dimers bound with both heads (D) , one head 
only (1) and empty sites (0) (see Fig. ||). We denote the 
probabilities to find a certain lattice site in one of these 
states by 2ud, n\ and no, respectively. Of course, nor- 
malization of the probabilities requires no+ni-\-2ri£i = 1. 
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FIG. 4: A configuration of bound kinesin molecules on a 
protofilament. The element D represents a dimer with two 
bound heads, the element 1 a dimer bound on one head and 
the element an emwty lattice site (vacancy). 
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A. Particle-hole symmetry 

Important information about the steady state and the 
dynamics of the system can be gained already by exploit- 
ing the symmetries of the kinetic process. Fig. || reveals 
a "particle-hole" symmetry by showing a transformation 
that maps the system onto one with the same reaction 
scheme, albeit with transformed kinetic constants. The 
details of the transformation are listed in Table 



the symmetry relation 

v(l/K 1 c,K 3 )=2-v(K 1 c ) K 2 ). (3) 

The symmetries can be best seen in a logarithmic-linear 
plot as shown in Fig. |(| From these relations we already 
conclude that the stoichiometry of dimers bound with 
both heads no reaches its maximum at K\c = 1 and 
that the total stoichiometry is v = 1 at that point. 
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FIG. 5: The one-dimensional dimer decoration model 
is invariant upon exchanging empty lattice sites (0) with 
single-bonded dimers (1) and the reaction rates according to 
(k +1 c,k f 2± ) «-» (fc_i,fc| ± ). 



The "particle-hole" exchange operation replaces 
dimers bound on a single head ("particles") by vacant 
sites ("holes") and vice versa. Dimers bound with both 
heads are kept invariant under the transformation. To 
obtain a system with an equivalent reaction scheme the 
transition rates have to be transformed as well. This 
is achieved by exchanging the attachment and detach- 
ment rate of the first head, /c+ic <-> k—\, and simultane- 
ously switching the forward and backward binding rates 
of the second head, fcf ± <-> k 2 ±. As a consequence of this 
symmetry operations the equilibrium constant K\c is re- 
placed by its reciprocal value while the binding constant 
of the second head K 2 is left unchanged. 

This symmetry has immediate consequences on the 
equilibrium stoichiometry. Since the coverage in the 
steady state is solely a function of the quantities K\c 
and K 2 this symmetry implies that the mean number of 
dimers attached with both heads (which are mapped onto 
themselves) ud{Kic, K 2 ) is invariant upon interchanging 
the attachment and detachment rates of the first head, 



n D (l/K 1 c,K 2 ) = n D (K 1 c,K 2 ) 



(2) 



values 1 19 
to Ki 



B. Experimental parameter values 

The kinetic constants for the binding of kincsin on 
microtubules have been determined by several groups 
0, ||, ||, Hi- The binding constants for the first head 
in the pre sence of ATP and low ionic strength have the 
~ k+i = 20^M~ 1 s~ 1 and k-i = 25 s -1 , leading 
0.8 ^M -1 . k-i is much smaller in the presence 
of AMP-PNP or in the absence of a nucleotide, about 
0.01 s- 1 jn|. The binding constant K 2 can be estimated 
indirectly as the ratio between the detachment rate of 
the monomeric and the dimeric kinesin and has values 
between 2.7 (in the presence of ATP) and 20 (without 
a nucleotide) Q. These values show that the model 
parameters depend strongly on the chemical conditions 
and the type of motor protein used in the experiment. 
A theoretical investigation allows for a systematic analy- 
sis and detailed classification of all the different regimes 
of binding kinetics in such a broad range of parameter 
values. 



III. EQUILIBRIUM STOICHIOMETRY 

In this section we review results obtained for the equi- 
librium stoichiometry of the dimer binding model |p.2| . 
In this section we will discuss the one-dimensional model 
where some short-range interaction between the dimers is 
taken into account; two-dimensional models are studied 
in Sec. |v} 



A. Analytical solution for the binding 
stoichiometry 

For the one-dimensional model the value of the mean 
occupation numbers in the steady state can quite eas- 
ily be determined upon using detailed balance and the 
fact that the dimers have only a hard-core interaction. 
Because we have non-interacting dimers the probabil- 
ity to find a certain sequence of 0's, l's and Vs {e.g. 
"0,1,D,D,1,D" , see Fig. ||) has to be invariant against 
permutations of these states. In such a random sequence 
the probability pi to find a particular state 0, 1 or D at 
a certain site i is given by 



Similarly the mean total number of bound heads per lat- 
tice site (binding stoichiometry), v — 2(n\ +no), obeys 



Pi = 



n + ni + n D 



(4) 
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TABLE I: The dynamics of the one-dimensional dimer decoration model is invariant under the following "particle-hole" trans- 
formation. 





state 


kinetic constants 


equilibrium con- 
stants 


Original 
Transformed 


1 D 

1 D 


fc+ic k-\ k^_ 2 k h +2 k{_ 2 

k \ ]^ C fc_|_2 ^^-2 ^ 2 2 


K lC K 2 
l/(K lC ) K 2 



Detailed balance requires that for each pair of possible 
configurations, their respective probabilities are in the 
same ratio as the transition rates between them. Hence 
the ratio of probabilities to find a sequence with a f or 
at a certain place is 



Pi 
Pa 



fc+ic 



K lC . 



Similarly, we get for transitions between D and 01: 
uf 



VP 
PoPi 



•J 



ft +2 



(5) 



(G) 



These two equations, together with the normalization 
condition, uniquely determine the values n , n-y and rip. 
The number of dimers bound with both heads is given 
by 



1 1 



n D = 



1 



2 2V(l + ^ic) 2 
and the number of dimers bound with one head 
K lC ( iK^c 



K lC +l \(l + K lC ) 



1 



(7) 



(8) 



The stoichiometry, i.e. the total number of heads per 
binding site v = 2(ni + no) then reads 



K lC - 1 / \K X K 2 C 



K lC + l \{l + K lC ) 



(9) 



The number of dimers bound with both heads per lattice 
site reaches its maximum n™ ax = (1 — l/V-^2 + l)/2 for 
K\c — 1. For an illustration of all these equations see 
Fig. §. 



B. Alternative derivation of the binding 
stoichiometry 

Here we present an alternative derivation of the bind- 
ing stoichiometry which has the additional benefit that 
it can also be used for two-dimensional lattices where 
dimers can bind in either direction. The idea is to re- 
late the equilibrium stoichiometry of the flexible dimer 
model to that of a stiff dimer model - a model where 
dimers can only bind as whole. Once we have specified 
the mapping one can determine the stoichiometry of the 
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FIG. 6: Binding stoichiometry v as a function of the solution 
concentration c and the constant K 2 . The thick curves show 
the total number of heads per lattice site v (^). The thin 
curves show the fraction of binding sites occupied with double- 
attached dimers 2nr> (^). 



general model from the stoichiometry of the stiff dimer 
model in one dimension which has been known for quite 



a long time 
the solution 



I already and represents a special case of 
3y McGhee and von Hippel p3[ 



1 



1 - 



1 



(10) 



Note, that in this notation the binding constant K con- 
tains the concentration c. In equilibrium, a mapping 
between the stiff and flexible dimer model is possible 
because of detailed balance and because transitions be- 
tween states and 1 are uncorrelated with the surround- 
ing configuration. Let's denote the equilibrium binding 
constant, the dimer and vacancy density in the stiff dimer 
model by K, no and no, respectively. Then, detailed bal- 
ance gives the following relation between these densities 
and the equilibrium binding constant, K — fiu/n^. In 
the same way detailed balance implies Ki = no / {tlqUi) 
and -fTic = ni/no for the general dimer model. The 
absence of correlations allows us to group together the 
states 1 and of the general model into a single state 
which we identify with the vacant state of the stiff dimer 
model. This gives the mapping as summarized in Table 
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TABLE II: Mapping between the equilibrium states of the 
stiff and the flexible dimer model. 





Stiff 




Flexible 


Densities 


fin, 




n-D, no + 11! 


Binding constants 


K 




K lC , K 2 


Equilibrium condition 


K = 


np 

"o 


K lC = Ki = 

n ' ' noil 



From this table we can infer that the number of fully 
bound dimers per lattice site in the flexible dimer model 
is given as 



ud(Kic, K 2 ) = ud(K) = no 



K 2 K lC 



Xl + K lC ) 2 

and the number of dimers bound with one head as 



(11) 



n\(Kic, K 2 ) 



l + Kic 



no(K) 



Kic 
l + K lC 



, 1/ K 2 K lC 

1 ~ 2 nD { OT^cf 



(12) 



In the one-dimensional case we can insert Eq. ( pf]| ) and 
reproduce Eqs. (Q) and (§h. Note that the mapping re- 
duces the number of parameters by one leading to a quite 
significant simplification. Though the stiff dimer model 
has not yet been solved analytically in two dimensions, 
the model with a single parameter can easily be solved 
by Monte-Carlo simulation (e.g., J24J). The entropy of a 
two-dimensional lattice fully covered with dimers is even 
known exactly [E5L E6|. 



IV. DYNAMICS 

In this section we will learn that the relaxation towards 
the equilibrium state described in the previous section is 
by no means a simple process characterized by a single 
time scale but actually shows quite interesting anoma- 
lous kinetics. We find different scenarios depending on 
the value of the equilibrium binding constants K 2 of the 
second head. If K 2 < 1 there are comparatively few 
double-bonded dimers in the steady state. Consequently 
cooperative effects and correlations are less important 
implying that the equilibration process is exponential 
with the time scale given by the rates fc+ic and fc_i. The 
dynamics changes drastically if the second head is very 
likely to bind, K 2 3> 1. Then, the equilibrium state will 
consist mainly of dimers bound with both heads aligned 
with twice the period of the lattice. Their positions are 
correlated at a length scale given by the average distance 
between defects in the periodic pattern (either vacancies 
or single-bonded dimers). As we will detail in this sec- 
tion these correlations lead to a dynamics which is slowed 
down drastically as compared to the kinetics of individual 
dimers. 



In the following we will mainly concentrate on the 
case K\c <C 1, to which we refer as the "stiff dimer 
model", reflecting the fact that dimers are only likely 
to be found in the state with both heads bound. In 
this limit we can discuss the essential features of our 
model, which include the anomalous relaxation kinetics. 
Because of the particle- hole symmetry (Sect. [II A ) which 
transforms K\ c — > 1/ (K\ c) , the results apply equally well 
for K\c 3> 1. The other interesting limiting case is the 
situation with K 2 3> 1 and K\c as 1. In this limit the 
dynamics is essentially similar to the case K\c -C 1, but 
more complicated in detail because we have to consider 
all three states, i.e. empty sites, dimers bound with one 
head and dimers bound with both heads. 



The stiff dimer limit 



It turns out that most interesting features of the model 
can already be discussed if we restrict ourselves to the 
limit 



K 2 » 1, 



#ic< 1 



(13) 



with K = K\K 2 c fixed at a value which allows the equi- 
librium state to consist mainly of dimers bound with both 
heads. The first condition of Eq. ( |l3| ) states that the sec- 
ond head is unlikely to be found in the unbound state if 
it has a place where it can bind. The second condition 
implies that isolated single-site vacancies are unlikely to 
be occupied with single-bonded dimers. Therefore, our 
model reduces to a dimer deposition-evaporation model 
where the dimers can only bind and unbind with both 
heads at the same time. We will refer to this model as 
"stiff dimer model" in the following. The vacancy con- 
centration no in the steady state, given as 1 — n\ — 2no 
in Eq. (Q) and (H) then simplifies to 



= 1/VT+4K . 



(14) 



The stiff dimer model describes dimers that bind at 
once to a pair of empty lattice sites and is characterized 
by the effective attachment rate of a whole dimer to an 
empty pair of sites, fc+c, an effective detachment rate 
and a rate with which a dimer can move diffusively with- 
out detaching, r^. Depending on whether the diffusion 
or the detachment rate is larger, we will distinguish be- 
tween the diffusive and the non-diffusive case. The stiff 
dimer model is similar to conventional dimer deposition 
models p7| , however with the difference that the attach- 
ment and detachment constants fc + and k_ also depend 
on whether the neighboring sites are occupied or not (but 
we will show that this dependence becomes irrelevant in 
the non-diffusive case). 
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1. Mapping between the general model and the stiff dimer 
model 

In the following we describe how the kinetics of the stiff 
dimer model follows from our general dimer model. In 
the limit we are considering, the detachment of a single- 
bound dimer is always much faster than the attachment 
of a new one (k-i ^> fc+ic, following from K\c <C 1) and 
the attachment of the second head much faster than its 
detachment {kll 2 ^> k^! 2 , following from K 2 ^> 1). We 
do not put any other limitations on the reaction rates for 
now. 



Case (i): 



Case (ii): 



Case (hi): 



k p + ,k p _ 



. • -•" • /' • , • -^E A A A A 



A: 



(16) 



Case (ii): Front neighbor occupied. There are two 
pathways that lead to binding next to an occupied 
site. The first head can either attach on the site 
next to the occupied one or one site further away, 
both at the rate k+±c. In the first case the prob- 
ability that the second head will attach before the 
first one detaches is k b +2 / {k-\ + k +2 ). In the second 
case the probability that the second head attaches 
in front of the first one is fc£ 2 /(&-i + fc^ 2 + k b +2 ). 
These two terms taken together give the attach- 
ment rate to the given pair of sites 



k+ic 



k 



f 

+2 



k-i 



^+2 



fc_l 



"+2 



Ub 
+2 , 



(17) 



If the process of detachment starts with detaching 
the front head (rate kt 2 ), the probability that the 
whole dimer will detach afterwards is + 
^+2 + ^+2)- ^ ^ starts with detaching the rear head 
(rate k_ 2 ), this probability is k-\J{k-\- 
gives the total detachment rate 



This 



Case (iv): 



k»k D _ 



(~'-^m~'--^L'" ~~^mf ""^m ~~*^m~ ' ~~"^M~ '"'^m "'^m ''^m ~^^m"^^L 



k M = k-i 



k b 



k, J _ 



k-i + k b 



k-i 



J,f 1 h.b 



(18) 



FIG. 7: Four cases for the binding of a new dimer (grey) 
in a gap between already bound ones (black). The effective 
binding rates are k, ' ' ' and the unbinding rates k_ ' ' ' . 



In a general treatment, the attachment and detach- 
ment rates also depend on the occupancy of the neighbor- 
ing sites. We therefore have to distinguish between four 
cases: both neighbors empty, front neighbor occupied, 
rear neighbor occupied and both neighbors occupied. 

• Case (i): Both neighbors empty. A dimer attaches 
with its first head at the rate k + \C. Afterwards, the 
second head can attach at the total rate k* +2 + k b +2 , 
while the first head can detach with rate fc_i. This 
leads to the total attachment rate per lattice site 



+ k b +2 ) 



K +2 



(15) 



• Case (Hi): Rear neighbor occupied. This case is 
analogous to the previous one, except that the in- 
dices / and b have to be interchanged. The attach- 
ment rate becomes 



klc- 



k+ic 



k b 

ft +2 



k-x + k^_ 2 



k b . 




and the detachment rate 



k-i 



■J 



k-i + k 



f 



(19) 



(20) 



Case (iv): Both neighbors occupied. In this case the 
dimer binds into a gap of exactly two sites. If the 
first head attaches to the front site in the gap (rate 
fc+ic), the probability for the second one to attach 



before the first one detaches is k b i _ 2 /(k- 



^+2)- 



Together with the analogous case in which the first 
head attaches to the rear site inside the vacancy, 
the total attachment rate becomes 



The process of detachment starts with the detach- 
ment of the front or rear head. The detachment 
gets completed if the other head detaches too (rate 
k-i) and fails if the detached head reattaches again 



(rate k 



f 

+2 



^+2)- 



The detachment rate then be- 



k + c = fc +1 c 



ft +2 



k-i 



ft +2 



(21) 



After the front head has detached (rate k_ 2 ), the 
probability that the whole dimer will follow is 
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fe_x/(fc-i+fci 2 )- Adding the pathway starting with 
the detachment of the rear head gives the detach- 
ment rate 



k D = fc_! 



k_ 



k b 
ft +2 



k — i ^+2 



(22) 



Of course, all these rates obey detailed balance, which 
states that 



K = K\Kic 



(23) 



Processes in which one head detaches on one side and 
subsequently attaches on the other side also lead to an 
explicit diffusion of attached dimers. A diffusive step for- 
wards occurs if the rear head detaches (rate k b _ 2 ) and 
reattaches on the front side (probability k +2 l '(fc-i + fc+ 2 + 
fc+ 2 )). The hopping rate (which is, of course, equal for 
backward steps) then reads 



uJ i.b 
^+2^+2 



K 2 (k. 



-1 + 



(24) 



So far we have shown that in the limit K2 1 our 
model can be mapped onto a simple dimer model. But, 
depending on parameters, this can happen in two qual- 
itatively different ways. First, if the detachment rate is 
smaller than the diffusion rate L < r,j, which is equiv- 
alent to the condition 



k-i < 



1 



(25) 



K +2 



K +2 



we get the diffusive limit. In the opposite limit, a dimer 
is unlikely to make a diffusive step before it detaches and 
we call that the non-diffusive limit. In the non-diffusive 
limit, the effective dimer attachment (fc+c) and detach- 
ment (fc_) rates become equal in all four cases, i.e. they 
become independent of the surroundings of a dimer. 

To summarize, the stiff dimer model in the non- 
diffusive limit, to which we will concentrate in the fol- 
lowing, describes dimers that can bind to any free pair 
of lattice sites with equal rate, namely k + c and can de- 
tach with a rate The dimer cannot move along the 
lattice and it always stays joined (cannot regroup its two 
parts with parts of other dimers). Although dimers can- 
not diffuse directly, in the case of high lattice coverage we 
can still introduce an effective diffusion which describes 
dimers near a gap detaching and reattaching on the other 
side of the gap. This simplified model still captures all 
the essential aspects of the non-equilibrium dynamics. 



2. Related dimer models 

Similar dimer adsorption models have been studied 
previously p7L pi p|, 51. Privman and Nielaba p7t] 



studied the effect of diffusion on the dimer deposition 
process (neglecting dimer desorption). There are several 
key differences to our model, the most important of which 
is that diffusion without detachment results in a 100% 
saturation coverage, whereas a model with detachment 
leads to a limiting coverage whose value depends on the 
binding constants of the first and second head (see Eq. 
This also has important implication on the dynamics as 
discussed below. For example, as a consequence of a fi- 
nite coverage the final approach to equilibrium is not a 
power law but exponential and there are, in addition, in- 
teresting temporal correlations in the fluctuations in the 
steady state. 

Stinchcombe and coworkers |^8|, |2^] studied the effect 
of detachment on the adsorption kinetics but allowed for 
regrouping of attached dimer molecules (two monomers 
that belonged to different dimers during attachment can 
form a dimer and detach together). While such pro- 
cesses are allowed for some types of inorganic dimers, 
they are certainly forbidden for dimer proteins like ki- 
nesin where the linkage between its two heads is virtu- 
ally unbreakable. If regrouping is allowed the steady state 
auto-correlation functions for the dimer density shows an 
interesting power-law decay oc t^ 1 / 2 p^, ^9|. This behav- 
ior can be directly linked to the gapless spectrum in an 
associated spin model . If regrouping is forbidden this 
power-law decay is lost (and becomes an exponential to 
leading order) due to the permanent linkage between the 
two heads of the dimer. Intuitively this may be under- 
stood as follows: only if regrouping of dimers is allowed 
are there locally jammed configurations (Neel-like states, 
"101 Of 0" , with alternating occupied and unoccupied sites 
in which neither attachment nor detachment of dimers is 
possible) in the final steady state which slow down the 
dynamics. The difference between the occupation num- 
bers on even and odd sites then represents a conserved 
quantity slowing down the dynamics. In our model, the 
conservation law is trivial since the difference between 
occupation numbers on even and odd sites always disap- 
pears. 



In addition, we will see in Sec. IV A 4 that the autocor- 



relation functions of the dimer and the vacancy occupa- 
tion number show strong differences in shape and typical 
times scales of relaxation. 



3. Numerical analysis of the non-diffusive stiff dimer model 

To study the kinetics we choose the initial condition 
as typically used in an experiment, namely an empty 
lattice. Figure || shows simulation data for the aver- 
age vacancy concentration as a function of time for a 
set of binding constants K = fc + c/fc_ = KyK^c. We 
find qualitatively very different approaches to the final 
steady state depending on the value of the binding con- 
stant K. For if <§; I, where the off-rate fc_ is much 
larger than the on-rate fc+c, there is no crowding on the 
lattice and the dimeric nature of the molecules does not 
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affect the approach to equilibrium. Hence one gets, like 
for monomers, an exponential behavior with a decay rate 
k-. In the opposite limit, K 1, there is a pronounced 
two-stage relaxation towards the steady state. The va- 
cancy concentration as a function of time reveals four 
regimes, an initial attachment phase, followed by an in- 
termediate plateau, then a power-law decay and finally 
an exponential approach towards equilibrium. 




0.001 



FIG. 8: Vacancy concentration in the non-diffusive stiff 
dirtier model as a function of time t for K — k+c/k- = 
f, ...,10000. Dashed lines show the steady state concentra- 
tion (^|) for a given K. The line with circles shows the short- 
time limit for K = 10000 (p^), and the dashed line the result 
of the reaction-diffusion model (Eq. po[ ). 

a. Initial binding At short time scales, t <C kZ 1 and 
t -C rZ 1 , when only deposition processes are frequent but 
detachment or diffusion processes are still very unlikely, 
the kinetics of the model belongs to the class of prob- 
lems frequently referred to as random sequential adsorp- 
tion (RSA) models (reviewed by Evans |3f| and Talbot 
et al. ||2j). RSA Models have among others been used 
to describe the chemisorption of inorganic dimers like O2 
on crystal surfaces, the binding of reagents to organic 
polymers or chemical reactions on the polymer chain. 

The deposition of stiff dimers has first been studied by 
Flory p3| ]. Later, it was pointed out by Page Q (see 
also a note by Downton [p5[ ) that the final distribution 
depends on whether dimers bind with both ends in par- 
allel (" standard modeF) or with one end first ("end- on 
modeF ) . A more recent article addressing this point and 
giving analytical solutions for both cases was published 
by Nord and Evans J3(|. In the "standard model" the 
dimers bind with equal rates to all free pairs of lattice 
sites. In the "end-on" the first end of a dimer binds with 
equal rates to any free site. The second end then binds to 
either side if there is space, or to the only available side if 
the other one is occupied. This has the consequence that 
the probability that a dimer binds to either end of a gap 
is f.5 times higher than that it binds to a certain pair 



of sites in the middle of the gap. This obviously leads 
to some clustering of dimers reduces the number of gaps 
that remain in the end. 

It is even possible to exactly calculate the time- 
dependence of the lattice coverage J37|. In the "stan- 
dard model" the vacancy density during the RSA phase 
follows (Si §7) 



no(t) = exp (-2 + 2e~ k + ct ) 



(26) 



The vacancy concentration locks at an intermediate 
plateau 



plateau 




0.1353 



(27) 



This Flory plateau represents a configuration in which 
all remaining vacancies are isolated, causing the system 
to be unable to accommodate for the deposition of addi- 
tional dimers. In the "end-on" model the plateau vacancy 
concentration is 



plateau 
>0 



/ 2^(Erf(V2) - Erf(l/v / 2)) - f 



0.1233 , 

(28) 



slightly smaller than in the non-diffusive limit. 

When applying the RSA results to our model, we have 
to distinguish between the diffusive and the non-diffusive 
case. In the non-diffusive case, the binding rate at all 
pairs of sites, regardless whether their neighboring sites 
are empty (|l5|), one of them ([Tt], [f9]) or both ( pf| ) are 
occupied, is equal. This corresponds to the "standard 
model" . The time-dependent gap density is given by ( prjj ) 
and its plateau value by (^7|). 

In the diffusive case, the binding rates on both ends of 
an interval, ([f7]) and ([l9]), add up to a rate three times 
as large as the binding rate in the middle of an interval 
(|f5|), while the attachment rate for a pair of sites with 
both neighbors occupied ( pf| ) is twice as large. Therefore, 
the RSA phase in the diffusive limit corresponds to that 
of the "end-on model" . The plateau gap density is then 
given by (p8|). In both cases, the final configuration is 
independent of the asymmetry in the binding rates. 

b. Power-law regime The secondary relaxation pro- 
cess towards the final equilibrium state is much slower 
than the initial random sequential adsorption process. 
Starting out of a jammed configuration in the Flory 
plateau the dynamics shows a broad time domain with a 
power-law oc f~ x l 2 instead of a simple exponential decay. 
Similar multi-stage relaxation processes have been ob- 
served in dimer models with diffusion but no detachment 
p7| , the key difference being that the detachment process 
implies that the steady state has a finite vacancy density 
and the final approach to the steady state remains not 
a power law but becomes exponential. There are also 
interesting similarities. In particular, in both models a 
large portion of the final approach to the steady state is 
mediated by the annihilation of vacancies. This behavior 
can be explained by introducing a particle representation 
in the following way (analog to the adsorption-diffusion 
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FIG. 9: Time evolution of a state, consisting of attachment 
(fast) and detachment (slow) events. The coarse-grained in- 
terpretation includes only long-living states (in boxes). The 
effective steps include diffusion, pair annihilation and pair cre- 
ation (not shownV-i 
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model p7|). We denote each vacancy on the lattice as 
a "particle" A, and each bound dimer as an inert state 
(00). The detachment of a dimer then corresponds to a 
pair creation process 00 — * A A, and the decoration pro- 
cess to pair annihilation A A — > 00. Since we consider the 
limit i(>l, states with two vacancies A on neighboring 
sites have a very short lifetime. We may thus introduce 
a coarse-grained model by eliminating these states (see 
Fig. |). 

Then processes like A00 — > AAA — > 00A result in 
an effective diffusion for particle A with a hopping rate 
''hop = ft— /2 (fc_ is the rate of the first transition and 
1/2 the branching probability for the second) and an ef- 
fective step width of two lattice sites. The hopping rate 
is different for the last step before two particles annihi- 
late (AOOA — > 00 AA), namely fc_/3. Pair annihilation, 
AA — > 00, occurs with a rate ft+c, or 2fc_/3 from the state 
A00A. Pair creation, on the other hand, occurs mainly 
through the process 0000 -> 00AA -> AAAA -> A00A; 
the corresponding rate is k-/3K per lattice site and 
hence largely suppressed with respect to the annihilation 
process as long as the particle concentration is far above 
its steady-state value. These processes are summarized 



instead of algebraic. 

Models like the simple reaction-diffusion model A + 
A — y show interesting non-equilibrium dynamics |3§| , 
|m [Io| . One can show |], [| that asymptotically the par- 
ticle density riA(t) decays algebraically, n>t(t) oc t~ x / 2 , 
which nicely explains the slow decay observed in simula- 
tion data (see Fig. ||). Note that a mean-field like rate 
equation 



d t n A {t) oc -n\(t) 



(29) 



would predict ri^(t) oc t . In our analysis we can even 
go beyond the asymptotic scaling analysis and try to 
compare with exact solutions of the model for a random 
initial distribution with density p by Krebs et al. 
They find (adapted to our situation with two-site hop- 
ping) 



n A (t) 



1 

2^ 



da 

li u ( u(\ 



V^(2^) (< latcau ) : 



— 16ur hop i 



plateau 
l 



) + (< 



pin iieatA ' 



(30) 



Its asymptotic limit (first determined by Torney and Mc- 
Connel j§) reads 



n A = (32^r hop t)- 1 / 2 



(31) 



are in excellent 
30, Minor devi- 



in Tab. Ill 



Other processes involve more particles and 
are of higher order in terms of a power series in K . 
They are negligible if the number of vacancies is small, 
Ua *C 1, and the binding constant is large, K 1. In 
summary, for K 3> 1, our dimer model can be mapped 
onto a one-particle reaction-diffusion model A + A — ► 0. 
Pair creation processes, — * A+A, are highly suppressed 
during the first stage of the relaxation process. The give 
a significant contribution to the dynamics close to the 
steady state where the relaxation becomes exponential with a diverging characteristic time scale Tk oc K. 



Note that it is independent of the initial particle concen- 
tration in the Flory plateau riQ latca u . 

Our Monte-Carlo data (see Fig. ||) 
agreement with the predictions of Eq. 
ations at times between the plateau and the power-law 
decay are due to the assumption of a random particle 
distribution underlying the derivation of Eq. |3^. In re- 
ality the state after initial binding (RSA) contains some 
correlations which however do not reach beyond a va- 
cancy (Markov-shielding) In other words, the sizes 
of occupied areas are not exponentially distributed, but 
they are uncorrelated to each other and the higher-order 
correlation functions decay as fast as 2-point correlations 
(super-exponentially) . Short ranged correlations in the 
initial configuration do have some effect on concentra- 
tion at intermediate times and this causes the deviation 
between the simulation data and the theoretical curve 
which, however, is not large. According to (4^, |^| cor- 
relations can affect the concentration in the asymptotic 
limit even if the pair correlation function is short ranged 
and the order-n correlation functions decay at a length- 
scale of 0{n), but this is not the case in our model. 

The results from the A + A — > model also become 
invalid for very long times where the particle concentra- 
tion comes close to its equilibrium value. In this limit 
the dynamics becomes scale- invariant Q. The particle 
concentration can be written in scaling form 



n A (t) = K- 1/2 h(k-t/K) 



(32) 
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TABLE III: The processes in the A + A «-» model onto which the stiff dimer model can be mapped in the limit K ^> 1. 

Process Rate Branching Total 
(1) prob. (2) rate 

(!) aaa ( 2 ) ', 1 kl 



J_ 

2 2 

2 2fc_ 

3 3 



Hopping A00 -4 AAA -4 A00 k 

AnulnlatKU) . lOOA Q ^A^A ^ | °° AA \ -f 0000 fc 2 2fc 

AA00 

(i) J 00AA I (2) (2) fc _ i fc_ 

AA00 | fc-+fc+ 3 3K 
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FIG. 10: Scaled equilibrium autocorrelation functions for 
C (i) = KC (t), C D (t) = C D (t) and C P (t) = KC P (t) , 
t — k-t/K. The function Cp(t) is compared with the an- 
alytical result (^), represented by the solid line. Simulation 
data were obtained with K = 100 (circles, for Co and Cd 
connected with dot-dashed and dashed lines) and K — 400 
(triangles) . j— i 
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4- Steady-state autocorrelation functions 

A limitation of the above mapping becomes evident if 
one considers the equilibrium autocorrelation functions. 
Contrary to conventional A + A <-> models there are 
three different autocorrelation functions with different 
functional forms (Fig. [H]). In the following we will use 
nui+x) as the occupation number (0 or 1) of a dimer on 
the pair of sites (i, i + 1) and ho i as the occupation num- 
ber of a vacancy at the site i (1 if there is a vacancy, 
otherwise) . The dimer autocorrelation function 



C d (t) = (n(i ji+ i)(t)n( M+1) (t + t)) - (h(i^ + i){t)Y 



(33) 



describes the probability to find a dimer at a given pair 
of sites simultaneously at times t and t + r. Similarly we 
define the vacancy autocorrelation function 



C (r) = (noMn 0i (t + r)) - (M*)} 5 



describing the probability that the site i is vacant at times 
t and t + r. To calculate their values at r = 0, we need 
the expectation values 

(n(i,i+i)(t)) - (n%, i+1) (t)) = n D = \{^- 



and 



{n Qi {t)) = (hl i {t)) = n Q 



1 



(35) 



(36) 



A third way of defining an autocorrelation function is to 
ask for the probability that a vacancy is either on the site 
i or on the site i + 1 

C p (t) = (^(h Qi + ho i+1 - n 0z n 0i+1 )(t) x 

x (h 0z +n Qi+1 - h Qi n 0z+1 )(t + r)\ 

- ((no, + h Qi+1 - n 0i n 0t+1 )(t)) 2 , (37) 

with 

(h Qi + ho i+1 - n 0i "n 



0i+i) - 7 



2n _ 2 
+ n ~ V4iTTT + 1 



(38) 



The important difference between Cp and Co originates 
from the fact that the vacancies diffuse by hopping over 
two lattice sites. This is why a vacancy diffusing along 
odd sites will not influence the function Co if the original 
vacancy was on an even site. The function Cp does not 
make this distinction. 

For K 3> 1 the autocorrelation functions become scale 
invariant as well. Their scaling form reads Co(ij = 
KC (t), C D {i) = C D {i) and C P (t) = KC P {t) with 
t = k-t/K. The latter corresponds to the autocorre- 
lation function in a reaction-diffusion model A + A «-> 0, 
which has recently been calculated analytically by Bares 
and Mobilia El 



C P (t) 



-21. 



(34) 



Erfc V 2t Erfc V 2t , 



with Erfc (x) = -= / e~ y dy . (39) 
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Note that the described solution is derived for a parame- 
ter set which requires a special relation between diffusion, 
pair-creation and annihilation rate to be fulfilled. In our 
model this is not the case. But, this difference becomes 
irrelevant in the scaling limit, since their and our model 
can be mapped onto each other by introducing a short- 
ranged interaction between particles. See also Ref. |4f| 
for a comment about the validity of the analytical ap- 
proximation used in Ref. plj ]. The observed deviations 
however do not affect the scaling limit. 

The other two functions decay on the same time-scale, 
but with different prefactors. The reason is that even 
if a pair of vacancies annihilates, the system still keeps 
memory on whether the surrounding dimers were located 
on even or odd locations and this gives those correlation 
functions that distinguish between even and odd sites a 
longer decay time. 




B. Two-particle model 

Another interesting case arises if we consider the limit 
K2 — * 00 with K\c being of the order of magnitude 1. 
In this case neither the dimers bound with one head nor 
the empty lattice sites dominate and we have to intro- 
duce two particle types. As we are used to, particles A 
should denote vacancies on the lattice. In addition, we 
introduce particles B, representing dimers bound with 
only one head. 

Reactions A — > B and B — > A occur at rates fci+c 
and ki—. If we assume that both attachment rates of the 
second head k£ + and k\ + are significantly higher than 
k\ + c and fci_, we are again dealing with the diffusive 
model ( pq ) . The diffusion mechanism for vacancies (par- 
ticles Ajis exactly the same as in the model with stiff 
dimers and according to (24) the hopping rate is given 
by Eq. (pi[). The diffusion rate of single-bonded dimers 
(particles B) has to be equal. This is due to the symme- 
try of our model upon exchanging particles A and B and 
the transition rates fci+c and k\^ (Tab. ||) and because 
the latter are irrelevant for the hopping rate. If a particle 
A and a particle B reach neighboring sites, they annihi- 
late quickly, while two particles of the same type do not 
interact. To summarize, we obtain a reaction diffusion 
model of the type: 



A- 



A^B 
B <-> 



While its dynamics is more complex at short times, the 
model becomes equivalent to the A + A <-> if the tran- 
sitions between A and B become faster than the typical 
annihilation time, which is given as the diffusion time 
between two particles, displaced by the average distance 
between particles on the lattice, ~ — Therefore the 

power-law behavior, given by ( |3l| ) remains untouched by 
the fact that we are dealing with two different particle 



FIG. 11: A computer simulation showing the time depen- 
dence of the vacancy (jia) density, the single-bonded dimer 
(n_s) density and their sum for K\c = 0.5 and K2 = 10 6 . The 
dotted line shows the predicted long-time limit of reaction- 
diffusion models given by Eq. (^) . In detail, parameters are: 



k+ic = 0.1, k-i = 0.2, kL = 10 6 , k f _. 



1, k 



+2 



0.1k 



f 



and kt 



0.1 k_ 2 . The simulation has been performed on a 



lattice with 2 sites and periodic boundary conditions. 



types. An example of a simulation in this regime is shown 
in Fig. 0. 



Dynamics of interacting dimers in one 
dimension 



In many biological systems the interaction between at- 
tached molecules plays an important role. For example, 
in the case of kinesin, there has been an observation of 
coexisting empty and decorated domains which can only 
be explained by an attractive interaction |l2] ] . Similar ob- 
servation has been made on actin decorated with double- 
headed myosin p6[ |. An interaction between the dimers 
can be introduced by assuming that a dimer is more likely 
to bind to a pair sites if one or two neighbors are already 
bound. The binding rate then becomes k+ = k+ = Ak+ 
(one neighbor bound) or k+ = A 2 k+ (both neighbors 
bound). Similar, we assume that a dimer with one bound 
neighbor dissociates with rate k^_ = k 1 ^ 1 = Bk^f and that 
a dimer with two bound neighbors dissociates with rate 
k D = B 2 k N . The cquilibr ium state of this model is still 
exactly solvable but the interaction changes both 
relaxation stages quantitatively. First, the vacancy con- 
centration on the intermediate plateau lowers since the 
interaction improves the formation of contiguous clus- 
ters during the first stage. The RSA phase can then be 
described in terms of Kolmogorov's grain-growth model 
p7|. The vacancy concentration after the initial bind- 

27^T @ E§' where r *»* is 



mg is given as no 
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the nucleation rate per free lattice site (in the simplest 
case, when K ^> 1, it is simply k + ) and r gr0 wth the sum 
of speed with which both boundaries of a nucleus spread 
over the lattice (normally k++k+ — 2Ak^). The second 
effect of the interaction is that the diffusional relaxation 
slows down by the factor B since the detachment rate 
decreases. And finally, the equilibrium vacancy concen- 
tration decreases |p3| . Nevertheless, interacting models 
show the same two-stage relaxation behavior. An ex- 
ample of a model with interaction is shown by the dot- 
dashed line in Fig. O. 



indirect interaction mediated through distortions of the 
underlying tubulin lattice or a combination of all of these 
possibilities. All of these mechanisms are consistent with 
the observation that the interaction is stronger on flat 
tubulin sheets than on cylindrical microtubules fjjfl. 
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FIG. 12: The number of vacant lattice sites as a function 
of time non-interacting (solid line) and interacting (dashed 
lines) dimers. The parameters were K = 100, A = 10 and 
B — 1/10 (a dimer is 10 times as likely to associate to and 10 
times less likely to dissociate from a certain site if one of the 
neighbors is present) 



V. DIMER MODEL WITH INTERACTION 



As we mentioned in Section IV C, interactions between 
bound dimers can often play an important role. In the 
kinesin-microtubule system, evidence for the existence 
of an attractive interaction comes from two key exper- 
imental observations. The first one is that bound kinesin 
dimers are found to form two-dimensional crystalline lat- 
Within the non-interacting model one can ex- 



tices 1 50 



plain longitudinal correlations along a single protofila- 
ment, but by no means the alignment of dimers on dif- 
ferent protofilaments. The second observation support- 
ing an attractive interaction between bound dimers is the 
coexistence of empty and decorated microtubules in one 
and the same experiment fl~2]| . Such a phase segregation 
is a clear indication that the system is in a coexistence 
regime with the strength of the dimer interaction above 
some critical value. The detailed nature of the interac- 
tion between the dimers is not yet known. It could be a 
direct interaction between kinesin heads and tails or some 



FIG. 13: Illustration of the interaction energy between two 
adsorbed dimers inrdnfferent relative positions. 
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Even if one assumes that there is only nearest-neighbor 
interaction, a general description would still require 12 
different coupling strengths (2, 4 and 6 for describ- 
ing the attraction between two single-bonded dimers, 
two double-bonded dimers and between a single- and 
a double-bonded dimer, respectively). Upon neglecting 
that the tubulin lattice is actually rhombic and not or- 
thogonal one can reduce the number of different coupling 
constants to 9. These are still too many for a general 
discussion. Hence, mainly for simplicity, we will restrict 
ourselves to a model with only two interaction constants, 
Jl acting in longitudinal and Jt acting in transverse di- 
rection, as illustrated in Fig. [l3|. We assume that the 
interaction between two double-bonded dimers on neigh- 
boring protofilaments is zero if they are displaced by one 
lattice site or more (Fig. [T^d). Otherwise, the transverse 
interaction is assumed to have the strength Jt between 
two double-bonded dimers (Fig. |l3]c) and ^ J t between 
a double- and a single-bonde d dim er (Fig. U3b). We use 
the same notation as in Sec. IV C| and introduce factors 
Ai and Bi measuring the change in the attachment and 
detachment rate by the presence of a neighboring dimer 
(with i indicating the relative position of the neighbors 
as shown in Fig. 13), 



7 with neighbour a t without neighbour 

K l,2+ ~ A i K l,2+ 

■ without neighbour 



v l,2- 



Bi k 



1,2- 



Detailed balance states that 



Ji = /csTln 



Aj_ 

Bi 



(40) 



For simplicity we further assume that the binding and the 
unbinding rate are always changed by the same factor, 
Ai = 1/Bf, this assumption is irrelevant in equilibrium, 
but it simplifies the kinetics. 

Fig. [li] shows the stoichiometry curves, equivalent to 
those in Fig. ^, for the interacting model. The most 
dramatic effect of the interaction is that it broadens the 
plateau where most dimers are double-bonded, giving a 
stoichiometry of one head per lattice site. If the interac- 
tion is strong enough (see e.g. the curve with J/ksT = 2 
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in Fig. |l4|), the stoichiometry curves show quite a steep 
rise from zero to the plateau at v = 1. This is quite in- 
dicative of a phase transition, albeit smeared somewhat 
out due to finite-size effects. 



v 1.0 




FIG. 14: Binding stoichiometry of the interacting dimer 
model (dashed lines) compared with the non-interacting (solid 
line). The binding constant of the second head is K2 = 10, 
the interaction strength J = Jt = Jl = fcsT (short-dashed 
line) and 2ksT (long-dashed line). 

To determine the critical point we make two further 
simplifications. First, we assume that w e are in the stiff 
dimer limit, K2 ^> 1, as we did in Sec. IV A. This stiff 



dimer model can be interpreted as a spin model with axial 
next nearest neighbor interaction, also known as ANNNI 
models, reviewed by Selke In terms of these ANNNI 
models, our model is a rather peculiar limit since it has an 
infinitely repulsive interaction between nearest neighbors 
and attractive interaction between next-nearest neigh- 
bors in longitudinal direction, as well as an attractive 
interaction between nearest neighbors in transverse di- 
rection. In the literature on ANNNI models mainly the 
opposite situation was studied since it leads to frustra- 
tions, while the "zero-temperature" state of our model is 
trivial, namely dimers aligned, let us say, the even sites 
on the protofilaments. This justifies an approximation 
which further simplifies the model, namely by assum- 
ing that dimers can only bind to even sites. Then the 
problem simplifies to a lattice-gas model with an effec- 
tive Hamiltonian 



n = 



E 



(—JLfl2ijn2i+2,j — JTn2i y jfl2i,j+l ~ M^2i,j) 



(41) 



where the coordinates (2i,j) run over all even lattice 
sites, fi2i,j are the occupancy numbers (0 or 1) on those 
sites and [i = ksThi(KiK2c) is the chemical potential of 
adsorbed dimers. The lattice gas model can be mapped 



onto the 2D Ising model |32|] with Hamiltonian 

W = (— J' L SijSi + i j — JtpSijSij+i — esij) (42) 



ho 



with exchange constants 



J't — -jJt , J'l — ~,Jl , 



and spin variables 



3 M 



= 2n 2l j - 1 . 



(43) 



(44) 



The spin s has the value +1 if it is parallel and —1 if it 
is antiparallel to the magnetic field e, 



/i + Jl + Jt 



(45) 



The critical temperature of the two dimensional Ising 
model was first determined by Kramers and Wannier [p3| . 
In the absence of an external field, the model has also 
been exactly solved by Onsager |p4| . The condition for 
the critical point reads 



sinh 



Jl 
2k B T 



x sinh 



Jt 
2k B T 



= 1 



(46) 



or J c = 1.76 k B T in the isotropic case {Jl — Jt)- 

The exact solution helps us to determine the critical 
point, i.e. the interaction strength needed to obtain a 
phase transition between an empty and a decorated state. 
The full form of the binding stoichiometry as a function 
of concentration, however, would be equivalent to cal- 
culating the magnetization of the two-dimensional Ising 
model in the presence of an external field, which has not 
yet been achieved analytically. Therefore, one has to 
rely on computer simulations to obtain the stoichiometry 
curves. An overview on Monte-Carlo simulations on lat- 
tice gas models can be found in Ref. |^|. The situation 
simplifies when one is far away from the critical point, 
i.e. when the coupling is much stronger than its critical 
value, J 3> J c . This corresponds to the low-temperature 
limit in the Ising analogue, where the spin polarization 
is saturated and points into the direction of the external 
field. The number of spins with orientation ] is given 
as n = NQ(e), where 9 is the Heaviside function. Back 
in the adsorption model, according to Eqs. ([l3|)-(fl"5|), the 
condition reads 



v = ®(ix + J L + J T ) 



(47) 



This expression can also be understood directly. The 
boundary of a decorated domain will normally consist 
of many edge elements (each one with three out of four 
neighboring sites occupied) and some corner elements 
(with two out of four neighboring sites occupied). If 
the corner elements are stable, the decorated domain will 
grow; if they are unstable, it will shrink. The stability 
condition of a corner element with one bound longitudi- 
nal neighbor and one bound transverse neighbor is given 
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by Eq. (j47|). Deviations from (|l7|) are possible due to 
finite size effects. This is because small patches have a 
higher curvature in the boundary and are therefore less 
stable. 



A. Dynamics of the interacting 2D model 



than the growth time, 



t 



growth ^ ^nuc 



(49) 



We have shown in Sec. IVA3t that the annealing of 
gaps between bound dimers can be slowed down sub- 
stantially as compared to the transition rates for single 
dimers. This effect is even stronger in a two dimensional 
system. Here linear domain boundaries emerge instead 
of point defects (gaps). These boundaries can annihilate 
when two of them join. But their diffusion is significantly 
slower than the diffusion of single gaps, roughly inverse 
proportional to the number of sites in the boundary. 

The nonequilibrium dynamics of decorating an initially 
empty two-dimensional lattice is an interesting and com- 
plex problem. It has first been studied by Becker and 
Doring |56| (reviewed in Ref. 57 ). A crucial concept is 
the critical nucleus, a certain number of adsorbed atoms 
(here dimers) necessary to form a stable nucleus able to 
grow and spread over the whole lattice. An analytical 
expression for the nucleation time can be given if the in- 
teraction is strong enough, and the critical cluster size 
becomes as small as four, three or two molecules. For 
larger critical cluster sizes a rough estimate says that 
the nucleation rate is proportional to the kinesin concen- 
tration to a power that is typically half the number of 
nucleus-forming units. 

There are different scenarios depending on the relative 
magnitude of the time to form a critical nucleus, t nuc , 
and its growth time, i gr0 wth, i-S. the time it takes such 
a nucleus to grow and cover the whole surface (see Fig. 
|l5|a) . j70| If the nucleation time is larger than the growth 
time, 



>t 



growth 



(48) 



the whole lattice will most probably be completely cov- 
ered with dimers before other nucleation sites can form 
anywhere else on the lattice. The final state is a defect 
free lattice generated from a single nucleus. This is con- 
sistent with experimental observations |50| . At the same 
time single-site nucleation also explains how empty and 
fully decorated microtubules can be found coexisting in 
one and the same experiment ||l2|| . When the decora- 
tion starts, nuclei will form on certain microtubules and 
totally decorate them, until the concentration of free ki- 
nesin in solution drops below the critical value, where 
the equilibrium phase is the empty lattice. A numerical 
study of the nucleation process Q gives us an estimate 
of the minimum interaction strength needed to fulfill the 
condition ( |48| ) as J w 3fcsT. This strongly indicates 
that there is a interaction between kinesin dimers which 
is of the order of 3 ksT. If the nucleation time is smaller 



there is a certain probability that a second critical nu- 
cleus forms on the lattice before the first one created 
had a chance to cover the whole lattice. The probabil- 
ity of such or higher order events grows with decreasing 
nucleation time. Then, the lattice will be covered by 
several domains, locally ordered areas which may be out 
of phase on the microtubule lattice. In order to finally 
reach a homogeneous decoration with only one predomi- 
nant domain a secondary process is needed which leads to 
a coarsening of the initial grain-like structure. This sec- 
ondary process includes domain wall wandering and an- 
nihilation induced by the detachment and re-attachment 
of kinesin dimers. As mentioned above, such annealing 
processes are expected to be much slower than the unhin- 
dered growth of individual domains. Fig. |l5| b summarizes 
the possible phases and dynamic regimes for interacting 
dimers decorating a two-dimensional lattice. 



VI. INTERACTING fc-MERS IN ONE 
DIMENSION 

Non-interacting fc-mers (fc > 3) are well described with 
the mean-field kinetics of the model kA — > which pre- 
dicts a gap density decay cx i^ 1 /^ -1 ) |6|, |58|. In the 
continuous limit, also called continuous sequential ad- 
sorption (CSA) or "car parking problem" the length of 
unoccupied space decays as cx 1/ log t in the mean-field 
regime §9|, |(J. 

The dynamics becomes more intriguing when a strong 
attraction between the fc-mers is introduced. An exam- 
ple of such a system is tropomyosin |H, a protein that 
plays a key role in the regulation of muscle contraction. 
In muscle each tropomyosin molecule binding on an actin 
filament occupies seven actin monomers (although other 
isoforms covering 5 or 6 monomers exist as well |3l] ) and 
thereby obstructs the action of myosin motors. The bind- 
ing of tropomyosin to actin is strongly cooperative - the 
binding constant next to a bound molecule is 1000 times 
larger than at an isolated place Jg^] . 

The initial binding phase is analog to the problem of 
interacting dimers and can again be described using a 
grain growth model |l4|, |3l], [l?]]. After the initial phase, 
gaps with sizes between 1 and k — 1 sites will remain. 
Their sizes are distributed randomly with equal proba- 
bilities of 1/ik — 1) for each gap size. What follows is 
again a process on a much longer time-scale in which k- 
mers at edges of the gaps can detach and re-attach. If 
they detach from one side of the gap and then re-attach 
at the other side, the gap moves by k sites in one direc- 
tion (see Fig. |l6|a). The rate of such diffusive steps is 
given as the detachment rate of a molecule on e.g. the 
minus side of a block multiplied by the probability that 
a molecule re-attaches on the plus side of the neighboring 
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(i) 



Multi-site nucleation 




(annealing: slow) 




10 F 



10" 



(ii) 



Single-site nucleation 




10" 



10 




Decorated 



(fast) 



J/k B T 



a) 



b) 



FIG. 15: a) Computer simulation of (i) multiple-site nucleation with subsequent annealing and (ii) single-site nucleation. 
Both processes can lead to a defect-free decoration. However, annealing can be extremely slow whereas single-site nucleation 
immediately leads to the ordered state. The parameters in the simulation were the following: (i): K\K<xc = 0.1, Al = At = 3, 
Bl = Bt = 0.3 (ii): K1K2C = 0.001, Al = At = 10, Bl = Bt = 0.1. b) The solid line shows the phase transition above which 
a decorated phase exists. The dashed line shows the border between homogeneous (defect-free) decoration reached immediately 
through nucleation and decoration with domain walls. The curves were obtained from a computer simulation on a lattice of 
14 x 100 sites (nota-hhat the defect-free range would be somewhat larger on a smaller lattice). 
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block before another /c-mer can re-attach to the position 
where the first one has detached: 



^"hop 



= k 



M 



= k 



1 



1 



(50) 



Two neighboring gaps (see Figs. |l6| b and c), each contain- 
ing between 1 and k — 1 sites with randomly distributed 
probabilities, can either join to a single gap (Fig. [l6|b) 
or, if their total size exactly fits one fc-mer, annihilate 
(Fig. |l6|c). If the original gap sizes are g± and 32, the 
joined gap has size (gi + g%) mod k. The annihilation 
will take place in k — 1 out of (k — l) 2 cases, therefore its 
probability is l/(fc — 1). Any other gap size will be the 
outcome in k — 2 cases, therefore having the probability 



(fc-2)/(fc-l) 2 . 

Again we can represent gaps as particles A hopping 
randomly along the lattice and joining or annihilating 
when two of them meet. The pair creation processes 
become relevant only if the particle concentration ap- 
proaches its equilibrium value and will be neglected for 
now. Then, we can map our model to a diffusion- 
coagulation-annihilation model consisting of following re- 
actions 



A + A — ► A with probability (k - 2)/(fc • 
A + A^0 with probability l/(fc - 1) . 



1) 



The fundamental difference between the interacting k- 
mer and the dimer model is that the fc-mer model not 
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FIG. 16: Examples of effective reactions that occur after the 
detachment and attachment of a heptamer: Diffusive step (a) ; 
Pair coagulation A + A — » A (b); Pair annihilation A + A — » 
(c). 
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only contains pair annihilation but also pair coagula- 
tion processes, A + A — > A. According to Ref. [p3[ all 
diffusion-annihilation models of the type 2A — + IA (or 
generally mA — > L4.) belong to the same universality 
class. In the asymptotic limit, the "particle" concen- 
tration can then be borrowed from the exactly solvable 
A + A -> model [§ 0. Adapted to our fc-mer model 
it reads 



n A (t) 



1 



2 - Z ^87rr ho pi 



(51) 



with Z = (k — 2) / (k — 1) and fhop = fc 2 niop; the second 
relation results from the fact that in each diffusive step 
a gap jumps over k sites. Hence, we finally obtain 



UA(t) = 



1 



fc 2 x /27rr h o P t 



(52) 



Interestingly, the asymptotic particle concentration is in- 
dependent of its initial value, i.e. the final gap concentra- 
tion at long times does not depend on the intermediate 
gap concentration n° A . Note that it is even independent 
of the solution concentration c. 

Of course, the mapping to the A + A — ► IA model is 
valid only as long as the concentration of particles A is 
well above its equilibrium concentration. Otherwise, par- 
ticle splitting events like A — > A+A also become relevant. 
If the concentration c is high enough, the t^ 1 / 2 regime 
will be followed by a regime where, due to these split- 
ting events, the average gap size will finally approach 
one site, which then are randomly distributed over the 
lattice. The number of gaps will decay according to the 
mean-field prediction oc t^ 1 ^ k ~ 1 \ as it does in the non- 
interacting fc-mer model 
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FIG. 17: The time-dependent gap concentration (no) and 
fraction of empty lattice sites (no) for an interacting pentamer 
(k = 5), as obtained from the simulation. The dot-dashed line 
shows the prediction of Eq. ( p^ ) and the dotted line the mean- 
field power law oc i -1 '^ -1 '. The thin solid and dashed line 
show the equilibrium values of na and no using the expression 
from Ref. [|23|]. The parameters are: fc+ c = 100, fc+c = k+ c = 
1000, kl = l,k1 = k™ = 0.1, fc+c = 10000 and k° = 0.01. 



equilibrium, processes of pair creation — > A + A become 
important as well. Similar to the dimer model, these 
processes leads to a final exponential relaxation towards 
equilibrium. The equilibrium gap size and distribution 
of the interacting fc-mer model are known exactly from a 
study by McGhee and von Hippel pSf] . 



VII. SUMMARY AND OUTLOOK 

In this contribution we have discussed the kinetics of 
some macromolecular assembly processes relevant for the 
formation of functional structures in cells. In particular, 
we were interested in the dynamics of ligand-substrate 
binding, where the substrate is a one- or two-dimensional 
lattice and the ligands are dimers or oligomers. As exam- 
ples we picked binding of dimcric kincsin on microtubules 
and tropomyosin on actin filaments. There are other re- 
lated systems where protein ligands bind to biological 
macromolecules, e.g. proteins to DNA or antibodies to 
viruses. 

Our theory allows us to study the effects of steric con- 
straints, binding rates and interaction between neighbor- 
ing proteins on the binding dynamics and binding stoi- 
chiometry. Our key results are as follows. The collective 
dynamics turns out to be not only much slower than the 
kinetics of single molecules but it also shows some in- 
teresting dynamic anomalies. Quite generally, we find 
that the binding kinetics goes through several qualita- 
tively different stages. At small times when deposition 
processes are frequent but detachment processes are still 
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unlikely, the kinetics can be described in terms of sequen- 
tial deposition models. Such models have been analyzed 
in various contexts years ago so that there are many ex- 
act results available. After this initial phase the system 
is left in a "metastable" or "locked" configuration, i.e. 
a configuration with many small vacant regions between 
blocks of bound molecules prohibiting binding of addi- 
tional ligands to the substrate. As a consequence one 
obtains an intermediate plateau regime for the binding 
stoichiometry. For escaping from this locked configura- 
tion vacant regions have to merge in order to accom- 
modate for the deposition of additional ligands. These 
are extremely slow processes mediated by a sequence of 
detachment and re-attachment processes. Upon identi- 
fying the gaps with particles A we have shown that the 
dynamics in this regime can be explained by mapping 
it onto reaction-diffusion models. For stiff dimers the 
corresponding model is A + A —> resulting in a gap 
density that decays with a power law oc t" 1 / 2 before it 
approaches the final equilibrium value. Taking into ac- 
count that actual biological dimers have some level of 
flexibility the model had to be generalized to include a 
second species of particles B, representing dimers bound 
with only one head. The corresponding reaction-diffusion 
model is then of the type A <-> B and A + B <-> 0. 
While its dynamics is more complex at short times, the 
model becomes equivalent to the A + A <-> if the tran- 
sitions between A and B become faster than the typi- 
cal annihilation time. If, on the other hand, one intro- 
duces a strong diffusion of adsorbed dimers the power- 
law decay is preceded by a mean-field regime where the 
gap density decays with cx t —1 . The two-stage relax- 
ation and the mapping to the A + A — > model remain 
valid even if there is an attractive interaction between 
adsorbed dimers on the one-dimensional lattice. In the 
two-dimensional model with interaction, the sequential 
deposition stage goes over into a nucleation process and 
the relaxation stage into a much slower annealing process. 
Yet even then the two-stage relaxation remains qualita- 
tively the same. 

The annealing process for interacting k-mers with k > 
3 is fundamentally different from the dimer model. One 
finds that non-interacting fc-mers are always well de- 
scribed with the mean-field kinetics of the model kA — > 
predicting a decay in the gap density cx i — VC 5- Intro- 
ducing a strong interaction between the fc-mers, as is the 
case for tropomyosin binding on actin filaments, we find 
a mapping to the A + A — ► I A model, valid for concen- 
trations well above the equilibrium value, which shows a 
particle concentration decaying cx t~ x / 2 . Closer to equi- 
librium other processes become important as well. First, 
particle splitting, A — > A + A leads to an intermedi- 
ate mean-field regime similar to that of non-interacting 
fc-mers. Later the events of pair creation, — * A + A 
become relevant and lead the system into the final expo- 
nential relaxation towards equilibrium. 

All of the above results concern the statics and dy- 
namics of ligands decorating periodic structures in the 



absence of ATP hydrolysis. Such systems are "passive" 
in the sense that the ligands bind and unbind from their 
respective substrates but do not show any motor activ- 
ity, i.e. move along these molecular tracks. As described 
above, this allowed us to give a detailed description of 
decoration experiments and to analyze corresponding ex- 
periments quantitatively. We also were able to identify 
and quantify cooperative effects between kinesin dimers. 
As our results show, the relaxation time in many ex- 
perimental or biological systems can be extremely long. 
Therefore, one has to be aware that an experiment where 
the measurement is taken soon after the start of the deco- 
ration often does not show the equilibrium configuration, 
but rather some intermediate state from the relaxation 
process. 

A natural extensions of this work is to investigate "ac- 
tive" systems composed of an ensemble of motor proteins 
and their respective substrates in the presence of trans- 
port along the molecular tracks driven by the chemical 
energy of ATP hydrolysis. Understanding such systems 
lies at the heart of many cellular processes. Although 
single processive motors like kinesin or myosin V can 
move loads over considerable distances, the number of 
motors acting simultaneously when transporting vesicles 
through the cytoplasm can be pretty large. For instance, 
in Ref. p5l the force acting on a bead in the cell has 
been estimated as 100 — 200 pN, which implies that at 
least 20 — 40 motor proteins were involved. How these 
active motors cooperate at high densities and what role is 
played by interactions between them, is to a large extend 
a still unexplored field with many open questions. 

An even broader range of questions arises when one 
tries to understand how cells organize their interior to 
fulfill their various duties. While organizing its interior, 
the cell has to physically separate molecules or molecu- 
lar aggregates from each other, define distinct functional 
regions, and actively transport molecules between these 
regions. Such processes rely on the assembly and ordering 
processes of macromolecules inside the cell and on forces 
that may be generated by molecular assembly or by the 
action of motor proteins. Understanding the principles 
underlying those complex cellular processes not only re- 
quires biological and biochemical studies but also studies 
of physical processes. In this respect microtubules and 
motors have been used as in vitro model systems to study 
various aspects of regulation and self-organization of cel- 
lular structures. It was shown theoretically and exper- 
imentally that the dynamic instability of microtubules 
in combination with microtubule polymerization forces 
is sufficient to provide a microtubule organizing center 
with a mechanism to position itself at the center of a 
confined geometry p6| , |67[ . In a simple in vitro system 
consisting solely of multi-headed constructs of the mo- 
tor protein kinesin and stabilized microtubules the for- 
mation of asters and spiral defects was shown |68[ . By 
varying the relative concentration of the components a 
variety of self-organized structures (patterns) where ob- 
served. The pattern formation observed in these simple 
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in vitro systems are only first examples of what will be a 
much broader range of self-tuned or regulated structural 
and temporal organization in biological systems. The 
control parameter determining the structure was the rel- 
ative concentration of the components. In cellular pro- 
cesses there are, however, other possible parameters for 
controlling and regulating such as binding of a host of 
associated binding proteins. We expect that these and 
other related biological processes are a source for fasci- 
nating cooperative effects and a multitude of interesting 
nonequilibrium dynamic phenomena. 
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